CMB POWER SPECTRUM ESTIMATION 
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We explore power spectrum estimation in the context of a Gaussian approximation 
to the likelihood function. Using the Saskatoon data, we estimate the power aver- 
aged through a set of ten filters designed to make the errors on the power estimates 
uncorrelated. We also present an improvement to using the window function, Wi, 
for calculating bandpower estimates. 



Estimates of parameters, 9 pi from data will in general have correlated 
'e p e p ') is not necessarily diagonal. A Taylor expansion of 
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the log of the likclhood function, In £, about the values of the parameter that 
maximize it, 9* identifies C p with the inverse of the second derivative of ln£. 
The expectation value of this second derivative is an important quantity known 
as the curvature matrix or Fisher matrix, F: 
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where Ct and C n are the theory and noise covariance matrices, (AA }_= 
Ct + C n , where A is the data (notation is more thoroughly explained inld). 

Since a "»„ p is approximately equal to its expectation value, the parameter 

covariance matrix, C p ' , is approximately the inverse of the Fisher matrix. 

Knowing C p is useful for two different purposes. Firstly, it is necessary if 
power spectrum estimation is to be used as a meana.of "radical data compres- 
sion" . We have tried this with the Saskatoon dataa and the parametrization 
Ci = q B C^\ where C x = 1(1 + 1)Q/(2tt) and Cf^ refers to standard, untilted 
cdm normalized to erg = 1 and restricted to t within band B. Having estimated 
qs for ten contiguous evenly spaced bands from I = 19 to £ = 499 we can ap- 
proximate the likelihood function of 9, where 9 is some other parametrization 
of the spectrum: 
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where qsiS) = (Ci(8^ b / (Cf dm ) b and the brackets mean a logarithmic average 
across band B. Seffl for how well this Gaussian approximation works when 
complicated slightly by a marginalization over calibration uncertainty. 

The second use of C p (or equivalently the Fisher matrix) is for the vi- 
sual presentation of the power spectrum. We can plot linear combinations of 
power averaged through a set of filters where the filters are designed to pro- 
duce uncorrelated estimates. One particularly useful set of filters comes from 
Cholesky decomposition], which is simply LU decomposition for a symmetric 
matrix; find L such that F — LL T . Notice now that L^ 1 does diagonalize F 
since L~ 1 F{L~ l ) T is equal to the identity matrix. See L in the left panel of the 
figure. The transformation affects the parameters by taking q to Q = L T <\. To 
convert the estimate of each Qp into a bandpower estimate in band /3, (Ci)/s 
(plotted in the right panel of the figure), divide it by the sum over the fil- 
ter function, fg = LB/3/(Cf dnl )B- To find the bandpower prediction of another 



theory, (fiDp, average it over the filter function: (C\)f3 = YIb /b v~*)b/ Sb /, 
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Note that fg is playing the role of Wi/l in the usual bandpower procednrj 
One can use F 1 / 2 instead of L which has been done with the COBE dataa. 
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Left panel: Cholesky decomposition of the Fisher matrix. For each value of /3 the values 
of Lb (3 at the £ value corresponding to the center of band B, are connected by straight 
lines. Right panel: Estimates of the power spectrum from the SK data set as given by 
the observing team (pentagons), and the quadratic estimator (triangles). The quadratic 

estimates are uncorrelated because they are estimates of power averaged over the filters J g . 
The power estimates in the highest three bands have been averaged together. 

We estimated q with a quadratic estimator, which can be derived from 
a Gaussian approximation to the likelihood function, which the central limit 
theorem tells us is good in the limit of large data sets. The Gaussian approxi- 



mation is equivalent to truncating the Taylor series expansion of ln£(# + 59) 
after the 2nd order term. Doing so allows us to solve for 89 that maximizes 
the likelihood: 
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Note that due to the matrix inversions this is an order of n 3 operation, where 
n is the number of pixelized data points. Approximations to the weights (Ct + 
C n )~ l are necessary to make this estimator practical for very large data sets. 

If we restrict ourselves to map-making experiments with no constraint 
removals, and make the parameters the Ci 's then this reduces to the quadratic 
estimator independently advocated by M. TegmarkQ. 

The dependence of the right hand side on 9 p suggests an iterative approach. 
The estimation of a% for standard cdm took only three iterations starting from 
erg = 1 and converging on as — 1.41±0.08 (c.f. 1.43±0.08 via the full likelihood 
analysis). Therefore, for the power estimates shown in the figure, we used the 
quadratic estimator with Ct for erg = 1.41 standard cdm. 

The bandpower expected from a given theory, Cj, for a given experiment 
with window function Wi is usually calculated by (^ z QWi/l)/ J2i Wi/la. The 
optimal (minimum variance) filter, however, is not Wi/l but instead Y],, Fw 
where the parameters of this Fisher matrix are Ci and it is evaluated with values 
of C\ consistent with the data and/or prior information. Heuristically, this is 
an inverse variance weighting. We recommend that observers reporting single 
bandpowers also report these filters which will improve the bandpower method 
as a means of "radical compression" . In the limit that Ct is proportional to 
the identity matrix the two filters are equivalent: ^ ; , Fu> oc Wi/l. 
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